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Abstract 



We develop a model of tax evasion based on the Ising model. We augment the 
model using an appropriate enforcement mechanism that may allow policy makers 
to curb tax evasion. 

With a certain probability tax evaders are subject to an audit. If they get caught 
they behave honestly for a certain number of periods. Simulating the model for a 
range of parameter combinations, we show that tax evasion may be controlled ef- 
fectively by using punishment as an enforcement mechanism. 
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1 Introduction 

Despite significant progress in expanding AUingham and Sandmo's (1972) seminal theoret- 
ical paper to compensate for its weakness in predicting sufficiently high levels of individual 
tax compliance, the way enforcement mechanisms work from a dynamic perspective still 
needs further attention. Specifically, a new line of research is pursuing a multi-agent 
based simulation approach (MABS), to analyse how different enforcement measures influ- 
ence the compliance behaviour of individuals in complex systems, where possibly different 
agents interact. A literature review of MABS-models regarding income tax evasion can 
be found in Bloomquist (2006). 

We find it especially worthwhile to investigate AUingham and Sandmo's theoretical 
prediction that both an increase in audit probability and an increase in tax penalty 
enhance tax compliance. We perform our analysis in the light of group influence on 
individual behaviour and flnd that both measures work to reduce tax evasion. Our study 
is related to two other MABS-studies: Davis, Hecht and Perkins (2003) (DH&P) and 
Bloomquist (2004), who develops a powerful tax compliance simulator (TCS), adressing 
various possible triggers of tax evasion. Like DH&P we assume that tax payers who 
are audited become honest upon audit. While Bloomquist assumes that an individual 
becomes more risk averse when seeing someone in her network being audited, we equip 
caught tax evaders with a memory. It reminds them to remain honest for a certain number 
of time periods after non-compliance was detected. We interpret the enforced period of 
honesty as a punishment. The severity of punishment obviously depends on the number 
of periods detected tax cheaters need to remain honest for. 

The importance we accord to group influence regarding an individual's decision whether 
to evade or not mainly stems from the work of DH&P, who stress that viewing others 
evade a previously honest tax payer becomes susceptible to evasion with some probability. 
The study of Korobow et al. (2007) also illustrates that especially group effects are im- 
portant for an individual's decision whether to evade or not. They flnd that the existence 
of social networks diminishes compliance. 

To further analyse how tax compliance develops over time when individual decision 
making is subject to group influence we use the Ising model, a simple model from physics 
describing how particles interact under different temperature levels. We flnd this modeling 
framework particularly appropriate, because it attaches a large probability to a state in 
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which an individual takes on the type that dominates her neigbourhood. To conform with 
the Ising model, we assume that only two types of individuals exist, honest citizens and 
tax evaders. 

2 The model 

We use a 1000 x 1000 grid square lattice, where in every time period each lattice site is 
inhabitated by an individual (=spin Si) who can either be an honest tax payer Si = +1 
or a cheater Si = —1 , trying to at least partially escape her tax duty. It is assumed 
that initially everybody is honest. Each period individuals can rethink their behaviour 
and have the opportunity to become the opposite type of agent they were in the previous 
period. 

The neighbourhood of every individual is composed of four people, agents to the north, 
west, east and south. Each agent's social network may either prefer tax evasion or reject 
it. Various degrees 

of homogeneity regarding either position are possible. An extremely homogenous 
group is entirely made up of honest people or of evaders. No majority regarding either 
position exists only when the neighbourhood is completely mixed up. This is the case 
when two individuals repectively prefer each position. 

Individual decision making depends on two factors. On the one hand, the type of 
network every agent is connected with exerts influence on what type of citizen she becomes 
in the respective period. On the other hand, peoples' decisions are partly autonomous, 
i.e. they are not influenced by the constitution of their vicinity. The autonomous part of 
individual 

decision making is responsible for the emergence of the tax evasion problem, because 
some initially honest tax payers decide to evade taxes and then exert influence on others 
to do so as well. How large the influence from the neighbourhood is can be controlled 
through the "social temperature" parameter, T (units: J/ks)- 

Total energy is given by the Hamiltonian H = — J2<i,j> JijSiSj — BJ2i Si- The sum 
runs over all nearest neighbour pairs of spins. Si and Sj. Jij is the coupling between 
spins, which we assume to be constant {Jij — J) for all neighbouring spins. B is a 
positive parameter. It denotes the importance of the magnetic field for 
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the total energy; in the tax evasion context it could be interpreted as the influence 
of mass media, but we do not use it here (i.e. B — 0). The magnetisation is given by 
summing the corresponding vahies of all spins (= J2i Si). 

le — Si - Sj^, multiplied by — J, denotes the energy resulting from the interaction 
between the considered individual and her four closest neighbours. It is calculated by 
adding up the products of the respective individual's type (spin Si) and the type of 
each of her four neighbours (spins Sj). It is known since decades that for T > — 
2/ln(l + \^) ^ 2.2 half of the spins are +1 and the other half —1, while for T < Tc there 
is a majority for one direction. 

We simulate the Ising model by performing a spin-flip only if a random number 
between and 1 is smaller than the normalised probability (p) of a spin-flip: p — 
exp{-AE/kBT)/{l + expi-AE/ksT). 

The following table illustrates the probabilities of a spin-flip, given 

a range of possible structures for the neighbourhood and the different temperature 
levels we used in our simulations (cp. figures 1 to 4). 





Temperature T 


T = 0.25 


r = 2.0 


T = 2.5 


T- 3.0 


r = 25 


le 


Probability of a spin-flip 




~ 1 


0.982014 


0.960835 


0.935031 


0.579325 




1 


0.880797 


0.832019 


0.791392 


0.539915 


h^o 


0.5 


0.5 


0.5 


0.5 


0.5 


Ie = 2 


^ 


0.119203 


0.1679815 


0.208608 


0.460085 


/e = 4 


fti 


0.017986 


0.0391655 


0.064969 


0.420676 



The higher the level of the temperature is, the 

more a probability near p — 1/2 is accorded to a spin-flip, regardless of whether it 
implies a reduction {Ig — —4 and Ig — —2) or an increase in energy (Ig — 4 and Ig — 2). 

On the other hand, the lower the temperature is (e.g. T = 0.25), the more certain it 
becomes that flips, which cause energy to fall, take place and 
that flips, which cause energy to rise, do not occur. 

Finally, in the intermediate case where energy remains unchanged by a spin-flip (/g = 
0), the probabihty of a flip is equal to p = 1/2, regardless of the temperature level. 
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Applied to tax evasion we can interpret the model as follows: Tax evaders have the 
greatest influence to turn honest citizens into tax evaders if they constitute a majority in 
the respective neighbourhood. 

If the majority evades, one is likely to also evade. On the other hand, if most people in 
the vicinity are honest, the respective individual is likely to become a tax payer if she was a 
tax evader before. For very low temperatures, for instance T = 0.25, the autonomous part 
of decision making almost completely disappears. Individuals then base their decisions 
solely on what most of their neighbours do. A rising temperature has the opposite effect. 
The individuals then decide more autonomously. 

We further introduce a probability of an efficient audit (pa). If tax evasion is detected, 
the individual must remain honest for a number of periods to be specified. We vary this 
amount of time and denote the period of time that the caught tax evaders are punished 
for by the variable One time unit is one sweep through the entire lattice. Appendix 
A contains a Fortran source code for our model. 

For background information on the Ising model and how to simulate it, see Stauffer 
et al. (1988). 

Already FoUmer (1974) applied the Ising model to economics. 

3 Dynamics of the model 

We simulate tax evasion dynamics for various temperature levels and differently severe 

punishments (k). When k equals zero, no punishment is present and the model describes 

the baseline Ising model. Various degrees of punishment are introduced for different 

temperatures, by setting k consecutively equal to 10 and 50 periods for all considered 

levels of the temperature. The probability of an audit is sequentially increased in steps of 

one percent, from to 100 percent. For a given probability of an audit the dynamics of 

tax evasion (measured as portion of the entire population) is depicted over the range of 

300 time steps. Our three-dimensional illustrations thus depict 101 single time series, one 

for each possible level of the probability of an audit. Additionally we use two-dimensional 

illustrations to better convey the form of the tax evasion dynamics, when the probability 

^By adding the number of tax evaders and honest citizens together in one time period, taking into 
account the positive or negative sign every agent is marked with, one can also easily find the magnetisation 
for each time period. 
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of an audit is either at a realistic level {pa = 0.05) or at a rather high level {pa = 0.9). 

— Figure 1 goes about here — 

In figure 1 we set T = 25. If the penalty is high enough (e.g. k = 50), tax evasion can be 
reduced to 0% in the short run, given that the probability of an audit is sufficiently high. 
In the case of a penalty duration of 50 periods and a probability of an audit of 90% (left 
panel in the second row of figure 1), within only a few periods each individual eventually 
is compelled to remain honest. This happens, because spins fiip relatively often at this 
temperature, which is far above the critical level. The peaks we observe in the level of the 
tax evasion, result from the fact that 90 percent of the initial large number of tax evaders 
gets caught and after k iterations simultaneously becomes free to decide whether to evade 
or not. Roughly half of them choose to become non-compliant again after being regiven 
the opportunity to evade. Moreover consecutive peaks in non-compliance diminish less 
over time and it takes longer until perfect compliance is established the further out on 
the time scale the evolution of tax evasion is considered at. When allowing for more time 
to pass, one can see that evasion eventually does not hit the mark of zero percent any 
more. After around 8000 time steps an equilibrium level of about 2% non-compliance 
is attained, because the number of agents who can freely decide which type to take on 
stabilises at a level consistent with this portion of tax evasion. 

If punishment is set equal to 10 periods (right panel in the same row) tax evasion only 
approaches zero percent and finally comes to rest at a level of 9%. 

Because the length of punishment now is too small to reach full compliance, the peaks 
are somewhat greater than at k = 50: Additionally to those individuals who are released 
after k periods from having to remain honest there now also exist other individuals who 
may be tax evaders as wellj^ 

In the two corresponding time series plots where k is equal to either 50 or 10 periods, 

but the probability of an audit is much lower {pa = 0.05), one can describe the dynamics 

similarly as above. Obviously the probability of an audit now is too low to reach full 

compliance, even if A; = 50 (cp. row three of figure 1). These two pictures illustrate well 

^Some of the tax evaders who were not punished, because the duration of the punishment is too short 
to estabhsh full compliance, also remain tax evaders in the subsequent periods. 
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that punishment is a suitable enforcement mechanism when the probabihty of an audit is 
set to a reahstic level. The more periods individuals are forced to remain honest in the 
case of detection, the lower the resulting equilibrium level of non-compliance apparently 
is. For A; = 50 it is equal to 21% and for k = 10 the equilibrium level of tax evasion 
amounts to 39%. 

— Figures 2 and 3 go about here — 

We also consider two other temperatures above T^, T = 3 and T = 2.5. The lower the 
temperature is, the slower adjustment towards the equilibrium in the baseline model oc- 
curs (i.e. 50% non-compliance) As individuals become tax evaders more slowly the 
tax evasion problem is less pronounced already from the beginning compared to higher 
temperatures. But, because spins flip less frequently at lower temperatures, the same 
enforcement mechanisms may work less efficiently in the short run than at higher tem- 
peratures. When considering either of the time series with = 50 and pa = 0.9 (at T = 3 
or T = 2.5) one clearly sees that evasion cannot be reduced to zero percent any more. 
On the other hand, if the temperature is at 25 everybody becomes an evader within only 
a few periods, so that the enforcement mechanism quickly entails the entire population, 
given that the probability of an audit is sufficiently high. Yet, for the considered low 
temperatures T = 3 and T = 2, the enforcement mechanism does not encompass every 
person any more, because it takes longer that all individuals once take on the type of a 
tax evader. 

Thus full compliance cannot be established any more in the short run. On the other 
hand, when looking at the long-run, one can see that the equilibrium levels of tax evasion 
are the lower the smaller social temperature is. 

— Figure 4 goes about here — 

Finally, for the sake of completeness, we also introduce a temperature level below 

(T = 2). At such low temperatures individuals seldomly decide to become non-compliant, 

■^The amount of time necessary for this equilibrium of the Ising model to be reached under different 
temperatures can be read off in either of the three-dimensional illustrations in the figures 1 to 3, when 
considering the evolution of the tax evasion over time, given that p is equal to 0%. 
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because their vicinity which is mostly comphant on average exerts strong influence on them 
to be honest as weU. Therefore the equihbrium levels of tax evasion are smaller than at 
higher temperatures. Obviously, also for this temperature a higher degree of enforcement 
works to reduce non-compliance more. 

4 Conclusion 

Considering that individuals are likely to be influenced in their decision to evade taxes 
by their immediate neighbours, we found that regardless of how strong group influence 
may be, enforcement always works to enhance tax compliance. Both, a higher proba- 
bility of an audit and a larger punishment work together to enhance tax compliance. 
To exhaust the model's explanatory power regarding how group influence affects overall 
compliance, it appears interesting to perform a similar analysis for different initialisations 
(e.g. everybody is dishonest in the beginning). Also, it seems interesting to consider if 
and under what circumstances the system may show large fluctuations in tax compliance 
and whether these effects can be controlled, with the aim of flxing compliance at a high 
level. 

We would hke to thank Martin Hohnisch and Sabine Pittnauer for their valuable 
comments and their time for fruitful discussions. 
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Appendix A: Fortran source code for a single time series 

1 .'Input data: L= 1000 ,k=50,T=25 ,mcstep =300, p =0.9 

2 program tax.evasion 

3 implicit none 

4 .'Declaration of parameters and variables 
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L=1000 

both=L*L 

Lmax=(L+2) *L 

Lpl=L+l 

L2pL=L*L+L 

L2pl=2*L+l 



nteger , parameter 
nteger , parameter 
nteger , parameter 

nteger , parameter 

9 integer , parameter 

10 integer , parameter 

11 real , parameter :: T=25.0 

12 integer , parameter :: mcstep=300 

13 integer , dimension (L* (L+2) ) :: is 

14 real , dimension ( —4:4) :: iex 

15 integer , dimension(L*(L+2) ) :: mem 

16 integer :: ie ,m,mc , hon , i , j , ev , k , ibm=l 

17 real : : ex ,p , ip 

18 .'Probability of an audit (p) is set to 90% 

19 p = 0.9 

20 ip = (2*p-1.0) *2147483648.0 

21 '.Number of periods tax evaders need to remain honest if audited (k) 

22 k=50 

23 '.Initialisation : Set everybody to honest and their memory to zero 

24 do m=l,Lmax 

25 mem(m)=0 

26 is(m)=l 

27 end do 

28 '.Spin— flip probabilities 

29 do ie = -4,4,2 

30 ex=exp(-ie *2.0/T) 

31 iex(ie)=(2.0*ex/(1.0 + ex) -1.0) *2147483648.0 

32 ibm=ibm*65539 

33 end do 

34 write ( * , * ) p , , 

35 '.Dynamics of tax evasion is simulated over mcstep+1 time steps 

36 do mc=l,mcstep 

37 '.Set counter of honest individuals to zero 

38 hon=0 

39 do i=Lpl,L2pL 

40 '.First periodic border constraint 

41 if (i .eq.L2pl) then 

42 do j=l,L 

43 is ( j+L2pL)=is ( j+L) 

44 end do 

45 end i f 

46 '.Audited tax evaders must remain honest for k periods 

47 if (mem( i ) . gt .0) then 

48 mcm( i )^iiem( i )— 1 

49 is(i)=l 

50 else 

51 ie=is ( i) *( is (i -l)+is ( i + l)+is (i-L)+is ( i+L) ) 
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52 ibm=ibm* 16807 

53 if (ibm . It . iex ( ie ) ) is(i)=— is(i) 

54 end i f 

55 ! Counting the number of honest citizens 

56 if (is(i).cq.l) hon=hon+l 

57 ibm=ibm* 16807 

58 ! Each audited tax payer obtains a memory 

59 if (ibm . It . ip . and . is ( i ) . eq. — 1) meini( i )=k 

60 end do 

61 ! Second periodic border constraint 

62 do j=l,L 

63 is ( j )=is ( j+both) 

64 end do 

65 ! The number of tax evaders 

66 ev=both— hon 

67 !For the time steps 1 to 300 print the same quantities as for mc=0 

68 write(*,*) p , mc , ev/ real ( both ) 

69 end do !End time step (mc) loop 

70 end program tax.evasion 
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Temperature: 25 






Figure 1: In the first two pictures, the x, y and z-axis denote the probabiHty of an audit Pa (x), the 
considered time period (y) and the corresponding tax evasion portion (z), respectively. 

These two pictures thus show the tax evasion dynamics when holding the social temperature constant 
at 25 and 

when controlling for the number of periods (fc) that a detected tax evader must remain honest for. If 
either k or pa is equal to zero we get the standard Ising model. When k and Pa are nonzero and positive, 
the augmented version of the Ising model, i.e. our tax evasion model, applies. 

The remaining four pictures visualise specific time series, by holding additionally the probability of 
an audit at a constant level. We leave the length of punishment unchanged at either 50 or 10 periods and 
depict the tax evasion dynamics for two different probabilites of an audit, i.e. Pa = 0.05 and Pa = 0.9. 
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Figure 2: Tax evasion dynamics when the 

temperature is held constant at 3. 
The same simulation design as in figure 1. 
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Figure 3: Tax evasion dynamics when the 

temperature is held constant at 2.5. 
The same simulation design as in figure 1. 



Temperature: 2.0 






Figure 4: Tax evasion dynamics when the 

temperature is held constant at 2. 
The same simulation design as in figure 1. 



